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INTRODUCTION 



Despite decades of improvement in observations and experi- 
mental results, the nature of DM is still under investigation. 
It is often claimed that the accuracy of the measurement of 
the Cosmic Microwave Background spectrum up to the sec- 
ond acoustic peak provides very strong evidence in favour 
of CDM scenarios. Ho wever other DM mode ls, e.g. Warm 
Dark Matter (WDM) JSchaeffer fc Silk 1988^ ■ or scenarios 
in whic h DM has non-negli gible interactions with photons 
(IDM) jBoehm et al. 2002ct) . also predict an identical CMB 
spectrum to that measured so far. 

In contrast with CDM, however, WDM particles do 
have a cut-off in the linear matter power spectrum (P(fc);), 
at a scale fixed by the warmon mass (the linear matter 
power spectrum of IDM also has a cut-off but it is much 
less drastic than in WDM and fixed by the ratio of the elas- 
tic scattering cross section of DM with photons or neutri- 
nos to its mass). Beca use of this cut-off, WD M particles 
lighter than ~ 1 ke V iNaravanan et al. 20001 ) or 10 keV 
jYoshida et al. 2003h . which experience damping on scales 



ABSTRACT 

Dark Matter (DM) is generally assumed to be massive, cold and collisionless from the 
structure formation point of view. Although this translates into a scale-free power- 
law matter power spectrum, a more correct statement is that DM indeed experiences 
coUisional damping, but on a scale which is supposed to be too small to be relevant 
for structure formation. The aim of this paper is to present a Cold (although "coUi- 
sional") Dark Matter candidate whose matter power spectrum is damped, discuss its 
consequences on structure formation and see whether it is distinguishable from stan- 
dard candidates. To achieve this purpose, we first calculate the coUisional damping 
and free-streaming scales of conventional (namely neutralinos) and non conventional 
DM candidates (say light particles heavier than ^ 1 MeV but lighter than 0(10) 
GeV). The latter are supposed to annihilate to get the correct relic density and can 
be considered as Cold Dark Matter (CDM) particles in the sense that they become 
non relativistic before their thermal decoupling epoch. Unlike neutralinos, however, 
the linear matter power spectrum of light DM candidates can be damped on scales 
of ~ IQ^Mq (due to their interactions). Since these scales are of cosmological interest 
for structure formation, we then perform a series of numerical simulations to obtain 
the corresponding non linear matter power spectra P(k)ni at the present epoch. We 
show that because of small scale regeneration, they all resemble each other at low 
redshifts, i.e. become very similar to a typical CDM matter power spectrum on all 
but the smallest scales. Therefore, even if lensing measurements (at redshift below 
unity) were to yield a P{k)ni consistent with CDM models, this would not constitute 
a sufficiently robust evidence in favour of the neutralino to rule out alternative DM 
candidates, even if their linear matter power spectrum is damped on relatively large 
scales. 

smaller than ~ 10® Af© or ~ IC'M© respectively, would be 
excluded by the WMAP d etection of reionization at large 
redshift iKoeut et al. 20o3) should it be confirmed (the ma- 
jor caveat being the ef ficiency of reionization sources see 
e.g. iKnebe et al. (200a) ). On the other hand, it might be 
difficult to discriminate among typical CDM, WDM or any 
kind of particles whose linear matter power spectrum is 
damped at cosmological scales by using non-linear matter 
power spectra only, because of the small-scale regeneration 
mechanism which tends to make the WDM and CDM P{k)ni 
very similar, at least at redshifts z < 1. Thus, the determina- 
tion of the matter power spectrum by lensing measurements 
(Mcllicr ct al. 2002) would probably fail to constrain the na- 
ture of DM, although there might exist other properties (like 
the structural parameters of DM haloes) on which lensing 
observations may set tighter constraints. 



Although a warmon with mass just above ^ 10 keV 
is still a viable possibility, CDM remains the most popular 
scenario, mostly because: 

• it is generally considered as the most simple and natural 
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solution. Indeed, it does not rely on any particle physics 
parameter, in contrast with the WDM scenario where results 
depend on the warmon mass, 

• the most p opular CDM candidate - namely neutrali- 
nos l|Favet IATTT i - is still allowed by cosmology and particle 
physics while WDM candidates proposed so far, e.g. grav- 
itinos or sterile neutrinos, appear much more constrained 
and very clo se to be excluded u nless on e invokes peculi ar 
mechanisms jBaltz fc Muravama 2005Jiansen et al. 20021) . 

Nevertheless, it is worth to keep in mind that neutralinos 
(or any other well-known CDM particle candidate) have not 
been discovered so far, so there is room for other possibilities. 

Coming back to structure formation, DM particles are 
assumed to be fairly massive in order to avoid large free- 
streaming effects that would wash out the DM primor- 
dial fluctuations llDavis et al. IQs'Jl. Fxirt hcrmore. to pre- 
vent the Silk damping mechanism Jsilk 19 67 Silk 196^) from 
operating, one requires that they should not have any 
electrom agnetic interactions. This has actually led some 
authors ijOunn et al. 197^eebles 1982h to propose, in the 
early eighties, that DM could be made of Weakly Interact- 
ing Massive particles (WIMPs), a statement which is now 
commonly accepted. Because of this "weakly interacting" 
property, one generally neglects DM interactions in numer- 
ical simulations. As a result Dark Matter, and WIMPs in 
particular, are said to be coUisionless regarding structure 
formation, their linear matter power spectra being generally 
taken as a "true" scale-invariant spectrum with no cut-off. 

Obviously, such a statement is a bit too drastic. 
If they are produced thermally, the DM particles must 
have non-negligible interactions (weak or not, provided 
they are not electromagnetic) which might have conse- 
quences on structure formation. Those interactions turn 
out to be also very useful as they are at the basis of di- 
rect and indirect experimental detection o f DM candidates 
jCRESST collaboration 200]|Barrau 2001^ . So it is worth to 
try to assess their detailed impact on structure formation, 
instead of neglecting them a priori. 

DM elastic scattering cross sections turn out to be im- 
portant for structure formation because they fix the epoch 
at which Dark Matter thermally decouples, i.e. they spec- 
ify the epoch at which the coUisional damping stops and the 
free-streaming starts. It is therefore of crucial importance to 
determine the strongest interactions a WIMP can have and 
what are their effects on structure formation. This eventu- 
ally provides a definition of the particles really permitted by 
structure formation models, thereby extending the notion of 
WIMPs to all particles having interactions that are not nec- 
essarily weak but small enough so as to preserve fluctuat ions 
on scales where we know they exist jBoehm et al. 200 iD . As 
a result, this opens new possibilities for Cold DM candi- 
dates, bearing in mind that their P{k)i would then differ 
from a "true" scale invariant spectrum. 

The aim of this paper is to answer the question as to 
whether or not structure formation can distinguish stan- 
dard CDM candidates like neutralinos from other more ex- 
otic alternatives. To this purpose, we first review the dif- 
ferent damping mechanisms that can affect DM linear mat- 
ter power spectra, with emphasis on the "mixed" damp- 



ing mechanism introduced bv lBoehm et al. (20011) . We then 
compute the damping scales of neutralino primordial fluctu- 
ations and compare them with those of lighter DM particles 
(having a mass 0(10 GeV) > rridm > O(GeV)). The latter 
indeed turn out to be viable Dark Matter candidates pro- 
vided their dominant annihilation cross section are S-wave 
suppressed so as to satisfy both the relic density criteria 
and the gamma ray/radio constraints. This condition can 
be achieved for example by introducing a "new" kind of 
interaction (not based on the exchange of Standard Model 
gauge bosons) that would have eluded experimental searches 
up to now (Bochm & Favct 2003). 

The light DM P{k)i being potentially damped on scales 
relevant to structure formation, we perform four WDM sim- 
ulations ^ with different scale cut-offs and one CDM numer- 
ical simulation for the sake of comparison. We present the 
evolution of the non-linear matter power spectra with red- 
shift for all these models in section 4 where we perform a 
detailed comparison of our WDM models with the typical 
standard CDM spectrum. 

In light of these numerical results, we conclude that even 
if lensing observations favoured a standard CDM P{k)ni, 
they would not exclude the existence of a cut-off in the linear 
matter power spectrum, unless they could probe the distri- 
bution of matter up to very small scales with great accuracy. 
Therefore, we emphasize that it is really difficult to rule out 
"exotic" DM models on the basis of the non linear matter 
power spectra alone. 



2 DAMPING LENGTHS 

The objects we know {e.g. galaxies, clusters of galaxies) 
are thought to form by gravitational collapse of primor- 
dial density fluctuations. However, in order to grow signifi- 
cantly, these fluctuations need to survive well-known damp- 
ing mechanisms that would wash them out. If the domi- 
nant species that constitutes these fluctuations was a light 
neutrino or any relativistic species, then the relevant damp- 
ing mechanism would be free-streaming. This free-streaming 
would generate a cut-off in the linear matter power spectrum 
at very large scales and eventually turn out to be in contra- 
diction with small-scale observations since all galaxies would 
have to be created out of pancake and/or filament fragmen- 
tation. A cut-off at large scale is also expected if the main 
component of these matter density ffuctuations is baryonic 
matter, because of their very large interactions with pho- 
tons; this is known as Silk damping. Since both these models 
appear in contradiction with observations, there is a need for 
another kind of matter, called DM, that would be massive 
enough to avoid significant free-streaming and without any 
"standard" electromagnetic interactions. 



^ Light DM P{k)i are in fact expected to be oscillating. Also one 
expects some power at very small scales (albeit much less than 
in CDM scenarios). Therefore, they may appear intermediate be- 
tween standard CDM and WDM linear matter power spectra. 
However, we use a WDM P(k)i to perform our simulations so as 
to make the point regarding the small-scales regeneration mech- 
anism. 
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In this paper, we argue that — quite similarly to WDM 
scenarios where the free-streaming length can be as large as 
10^ — 10® Mq — there exist realistic scenarios where the col- 
lisional damping lengths of CDM can be significant (at least 
up to 10^ Mq). These DM particles would be Cold in the 
sense that they annihilate and decouple after they become 
non-relativistic, but on the other hand they would have a 
linear matter power spectrum very different from what is 
generally expected. We also show that another damping ef- 
fect, called mixed damping, can become particularly impor- 
tant in the case of light DM particles and is responsible for 
this peculiar behaviour. 

The latter can be understood as an intermediate ef- 
fect between collisional damping and free-streaming. It takes 
place when DM is thermally coupled to a species i that is 
already free-streaming. This means that it basically starts 
when species i thermally decouples (time tdec(i)) and lasts 
until DM stops being in thermal equilibrium with i (which 
occurs at time tciec{dm-i))- The condition for the mixed 
damping to take place is therefore written as: tdec{dm-i) > 

In the following, we shall focus on i = v, since DM- 
neutrino interactions are likely to satisfy the previous rela- 
tionship, i.e. 



^dec(^dm — u) 



> t. 



d&ci^f) 



or, equivalently. 



dec{dm — u) 



dec{u) 5 



where T'dec(dm-i') is the temperature at which the dm-v in- 
teractions stop having an effect on the DM fluid and Taec(v) 
is the neutrino temperature at their decoupling (say in a 
standard scenario: Tdec{,^) = T'tiec(i^-e) ~ IMeV). 

The mixed damping effect has never been considered 
before because its existence has only been recently pointed 
out bv [Boehrn_^t_aL_(2 00lfl (and its physical relevance by 
iBoehm fc Favet ('2003^ ) but also because DM candidates are 
generally thought to have thermally decoupled before or at 
the onset of 1 MeV {i.e. Tdm-v > Tdec(i/-e)) so that the 
condition for this damping to occur was thought to never be 
satisfied. 

Let us now suppose that we are in a scenario where 
Tdec{dm-i') < Taec{u)- The damping scale of DM primordial 
fluctuations then originates from i) the collisional damping 
effect due to the coupling of DM particles with species that 
are still collisional (this actually includes the neutrino con- 
tribution for T > TdEc(v)), ii) the mixed damping effect, 
originating from interaction with neutrinos and acquired 
during the period [ tdec(i'), idec(dm-^)l, and finally iii) the 
DM free-streaming for DM temperatures below rrfec(cim) (in 
fact Tciec(dm-v) if dm — v are the strongest DM interac- 
tions). The associated damping scales will be denoted /cd, 
Im d, Ifs respectively. Their expressions are given (according 
to lBoehm et al. ('200]|Boehm et al. f2002bl) ') by^: 



^ The normalization factors are not shown in these formulae but 
we do take them into account in our calculations. 



Pdm ^dm 

pi' dm. a? 



dt 



E 



.,2 

■ dt 



■ dt 



pi Vj 



*d<,c(,.) 



pHa'^ 



Ifs (X I — dt^Maxi — 

It a \ a I , 



(1) 

(2) 
(3) 

'dce(dm) ^'O 

with i denoting all the species that can transmit their colli- 
sional damping to DM fluctuations, Fdm the total DM inter- 
action rate, Fi = I'i-i the sum of the partial interaction 
rates of species i {j being the species in thermal equilibrium 
with i), H the Hubble constant and a = a{t) the scale-factor 
at a time t. Here tx denotes either ttiec{dm-i^) if idec(dm-v) < 
tdec(v) (in which case there is no mixed damping since DM 
thermally decouples from neutrinos before they enter their 
free-streaming regime) ortaec(,.) if idec(dm-,.) > idec(^) (since 
in this case the integration of l/F^ for t > tdec(i/) does not 
have any meaning). We emphasize that eq.lQ and eq.((2l 
supposes implicitly that taec{dm-i) < tdec(i) for i =^ or in 
other words, that the mixed damping species is only due to 
neutrinos. 

We also draw attention to the fact that we wrote Imd 
and led as two separate contributions when in fact, the mixed 
damping (eql^J is part of the collisional damping (eq0. 



Thus, even if t 



dec{dm — h') 



> tdec{v), we should write only one 



contribution, namely: 



1 

''cd{dm — f) ^ 

with 



pF^ a2 



dt 



pv c 
pHa^ 



dt 



~ / pTdm ^ ^cd(dm-x), 

but we will avoid this notation in what follows (and come 
back to eqQ and eq|2l to underline the importance of the 
mixed damping regime. 

The reason why neutrinos are expected to be the only 
species of interest for the collisional and mixed damping can 
be easily understood. Since the transport coefficients that 

2 

enter eq^are proportional to ''p"' , only photons and neu- 
trinos are expected to give a large contribution to led (and 
therefore to Imd)- Their velocity is indeed maximal and they 
do not experience any annihilations that could reduce their 
energy density. To contribute to the collisional damping, 
however, the DM-photon interactions should be such that 
tdec(dm--i) ^ tree (with tree the recombiuation epoch). Such 
a condition actually requires very large (and therefore un- 
likely) values of the dm — 7 elastic scattering cross section so 
neutrinos should eventually turn out to be the only species 
of interest as far as the estimate of led and Imd is concerned 
(since tdee{dm-v) ~ idec(v) is expectcd to rely on interactions 
of weak intensity). 

In what follows, we shall therefore make the reasonable 
assumption that tdec(dm--y) ^ tree, and neglect the photon 
contribution to the collisional damping effect. 
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Let us now compute the terms associated with neutrinos 
in led and Imd- The interaction rate — T^-j can be 
decomposed as = V^-e + F^-dm + J2j^e ^i^-i- We shaU 
neglect the terms '}2jj^e^'^-i ' since the Ti^-j are seen to 
be always smaller than T^-s, and make the assumption 
Ti^-dm < Fy_e. The neutrino decoupling is then given by 
Fy_e ~ H, as expected in a standard scenario. The non 
conventional case F^_dm. > Fi/_e will be discussed later for 
completeness. 



• Let us first study the the case where tdec(dm-v) < 



dec(iy) 



dec{u — e) • 



This is a standard situation without mixed damping ef- 
fects so we only need to discuss the neutrino contribution 
to led- Despite the term (? which appears very promis- 
ing, the maximum coUisional damping that standard neu- 
trinos can really transmit to DM fluctuations is l?cd{dm-v) J: 
j«d.c(.-=) _£M_i_^^i ^ [0(100 pc)]2 , (the maximum damp- 
ing length being reached when tdec(dm,-u) - td^c(u-^))- 
Thus, if a DM candidate decouples from neutrinos at ~ 1 
MeV, the cut-off in the linear matter power spectrum due to 
DM-neutrino collisions should be around O(O.lA/0), which 
is a remarkable property! 

This is actu ally much larger than the scale computed by 
iMisner (19671) (say 10"" - 10"^ Mq) relative to the damp- 
ing of electronic primordial fluctuations^. But one can check 
that his computati on, now rewritt en in terms of the Wein- 
berg formulation llWeinberg 197 j) . indeed corresponds to 



iMisner CX J ^^^(^ 

100 pc. 



dt , so that In 



This damping has to be compared with the free-streaming 
motion acquired by the DM during the period [tdec{dm), ia] 
{to denoting the present epoch and tdec(dm) = tdec{dm-v) ~ 
tdec(v-e) if dfTi — I' is indeed the last DM interaction). For 
annihilating particles heavier than a few MeV, which is a 
reasonable assumption in the case of a usual WIMP candi- 
date" decoupling at tdec(i^-e), we flnd a free-streaming mass 
of 

-3/2 



M 



fs 



rUdm 10"" 
MeV 10-1 



10^ Mq 



mdm 

MeV 



-3/2 



Mp 



SO Med ~ .1 f^)^^^ Mts fo r WIMPs decoupling as late 



as t 



dec{i/ — e) 



iBoehm et al. f200j) . 



• Let us now consider the case tdec(dn 

tdec(iy — e) • 

Since the mixed damping becomes relevant, the DM pri- 
mordial fluctuations are expected to be washed out by neu- 
trino free-streaming motion (instead of neutrino collisions). 
Free-streaming being extremely efficient, one expects a sig- 
nificant damping effect for a temperature Tdec(dm-v) just 

^ The existence of DM was actually not considered at that epoch. 
^ The damping scale of annihilating WIMPs, having a mass above 

a few MeV, is 1;, ^ 0(200) kpc (^)-''" {^^^^)^'\ 

or equivalently Af^, ^ W Mq (i^f)"'/' (^^^i^)''" , 
'^dec(dm) being the scale-factor when the DM thermally decou- 
ples. 



slightly below Tdec(u-e) (assuming standard neutrinos i.e. 
Tdec{u) = Tdec(u-e))- But cau Tdec{dru-u) be lower than 
Tdeciu) = Tdec(v-e) for Ordinary WIMPs? 



What determines whether the mixed damping regime 
exists or not is basically the strength of the dm — v interac- 
tions. The stronger they are, the lowest Tdec{dm-v) is, and 
the more likely the condition Tdec{dm.-u) < Tdec(y) is go- 
ing to be fulfilled. Accordingly, we now proceed to establish 
the relationship between Tdec{drn-v) and the interaction rate 



2.1 Relationship between Tdecidm) and Tdm-v 

The temperature can be written as Tdec{dm,-v) — 

Ta/adecidm-u) where ffldec(dm-i/) = {^dec{dm-v)/ Hr) 



is the scale-factor at tdec(dm-i^), ^dec(dm-i^) = 
Tdec{dm-^) a\ec(dm--y) is the effective (comoving) in- 
teraction rate at the same epoch and Hr ~ 2 10"^'^ s~^ (in 
the radiation dominated era, which is the relevant epoch 
for our purpose). 

The relationship between the interaction rate T(^dm-v) 
and the elastic scattering cross section oi^dm-v) between DM 
and neutrinos can be obtained by writing the Euler equa- 
tions fo r the DM and neu trino fluids (in the Newtonian 
gauge) iBoehm et al. 2002d) : 

vim = - Hvdm. - a fl {Vdm. - Vu) 

— '^k'^v — CI k (v 
-a /i {Vv — Vdm), 



with 

A = 0"(dm-i/) CUdm, k — (J{^^e) C , S ^ CC . 

Pdm 

Here H is the Hubble parameter, $ the Bardeen potential, 
5v the neutrino density perturbations, tt^ the anisotropic 
stress tensor and Vdm, Vv the velocity divergence of DM 
and neutrinos respectively. We then get: 

fdec(dm-i') = S ^ fl 'J"^!t^ec(dm-,x) (^) 

which implies, if we parameterize the cross section as 

'^(dm-v)C = b^K^dm-v) 

(n > ) and use the definition Tdec(dm-u) oc adec(dm-i.)- 

fdec(dm — t') ~ b^K^dm — u) '^i^'^ dec(dm—u) I '^dm ^ 

l/(n + 2) 



dec{dm — v) 

C^dec{dm — v) 
'^{dm-u)C 



To Hr mdn 

bel(dm — iy) ^ 

h ^ n + 1 \ l/(" + 2) 

OeUdm-u) n„ lO \ 

Hr mdm. ) 

h T" 

^el(dm — u) dec{dm~Lj) 



Kl(dm-u) ToHrmd^ 



(5) 
(6) 

(7) 



Here fii, = a is the comoving number density and 
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T/mdm a factor that reflects the efiiciency of the dm—v col- 
hsions in transmitting the neutrino momentum to the DM. 
In other words, our effective V(^cLm-v) corresponds to a stan- 
dard interaction rate (usually defined as the product of av 
times a number density) times the transferred momentum 
and T ciec{dm-v) is the same interaction rate, but comoving 
and estimated at tdeddm-v)- 

As expected, one can check from eq|3 that the 
largest the cross section (and incidentally &e!(tim-i/)) is, 
the lowest Tdec(,dm,-y) is- For n = 0, adec{dm-i^) ~ 

2.3 10*^ y/c^[dm~v)C (^v)"^'^^ (assuming that a^dm-v) in- 
deed yields a decoupling in the radiation dominated era). 



For n = 2 (the most likely value), we find: 

Tdec(dm-.) ~ 2.5 10 b^i' [^^) CI 



^ dec{dm — v) 



MeV. 

3 ^1/4 /nW\-l/4 
1/2 



(8) 



cr(dm-.)C -6 10 fej (,Mivj (9) 
(with bd = (&ei(dm-^)/ cm^ s"^)) 

2.2 Conditions for mixed damping to exist 

By using the definition^] we actually introduce a "non stan- 



dard" notation where, in fact, F, 



dec(i' — dm) 



7^ rdec( 



dm — h') • 



The relationship between these two interaction rates can 
be written as: rdec(^-dm) = ^j^'^ dec(dm-v) with pdm and 
Pv being the energy densities of DM and neutrinos respec- 
tively. Roughly, one can write Fj^.^m) = <^dm-v cndm and 

r(dm — i^) ~ ^dm — u Cfiu . 

If one requires neutrinos to decouple at ~ 1 MeV (which 
corresponds to the standard scenario) and DM to be still 
influenced by neutrinos, then the condition for the mixed 
damping effect to take place should in fact be rewritten (still 
assuming monotonic behaviour) as: 



- u — dm\t 
Fi^-elt 



£iec(i.-c) 



< 
< 



■ dm~h'\t , 



which means, with our notation: 



(10) 
(11) 



f"dn 



0"dn 



< 



ndr, 

"■d 

n 



T 



or (if one assumes that DM has finished annihilating before 
(or at the onset of) tdec{v-e))'- 



Cdn 



0"dn 



< 
> 



A cr^_ 



mdm Qe 



(12) 
(13) 



with A 



MeV 

^(27r)^/2 (C(3) the Riemann zeta funo 



tion of 3) and Xd = mdm/Tfo G [12 - 20] for mdm £ [1 - lO'^] 
MeV {Tfo being the freeze-out temperature). 

We will use the previous equations only for particles 



heavier than mdm ^ MeV (so that they stop annihilating be- 
fore primordial nucleosynthesis) . The mixed damping regime 
is expected to take place for DM particles with an elastic 
scattering cross section 



""(dm-i^) 

or 



> 2 10' 



r mdm \ 

\ MeV/ 



I, ^ in-56 / nidm \ 5-1 

beiidm-.) > 10 [-^) cm s 



(bearing in mind that (T(dm_^) cannot be too large so as to 
satisfy ea ll21 . 

The elastic scattering cross section of heavy DM parti- 
cles (with mdm ^ me) with neutrinos through the exchange 
of heavy particles X (with a mass mx ^ mw) is however 
likely to be such that (Jdm-v ^ o"y_e. Since this is in con- 
tradiction with eg 1131 one expects the mixed damping effect 
to be of interest only in non standard circumstances like 
resonant effects, light DM particles etc. 



2.3 Expressions of the damping lengths 

We can now express the different damping lengths in terms 
of the DM cross sections. Note that we assume n = 2, 
a-nr < adec(dm-v) < Heg, F^_dm < T^-e and we make the 
assumption that DM annihilates. 



cdi^dm — f) 



^dec{dm—u) \ / ^dm — u 



ffldec(i'-e) / yF(^_e) 



1/2 



\ 5/2 
^dec{dm — f) \ 



< 100 pc (10-'=^ b,j) 



MeV 



, J 1 ^dec(dm — iy) \ -r \ 

^md f'U — e I I it C['dec(dm — ^dec{u — e) 

\ CLdec(iy-e) J 

~ 100 PC (1.4 10^^ beif^^ ( '"''™ '^'^ 



hs(dm) ~ 200 pc 



Isd(dm) < 200 pc 



MeV 



7Tldm\~^/2 f (ldec{dm) 

MeV/ \adec(,.-e) 

mdm \ / Cldec(dm) 

MeV/ \adec{^-e) 



1/2 



1/2 



As a comparison, in a non-standard scenario where 

Tdec(v) = Tdec(v-dm) lS 1 McV Slid Tdec(dm-i/) < Tdec(i^-dm) ) 

the mixed damping would start at Tdf,c{i^-dm) instead of 
Tdec(iy-e)- This is important as most of the values of a(^dm-u) 
one could be tempted to consider may yield the scenario 
Tdec(^-dm) < Tdec(v-e)- ludccd, ouly a smaU range of val- 
ues of a^dm-iy) can satisfy the conditions ea llOl and ea llll 
associated with the mixed damping regime due to DM col- 
lisions with "standard" neutrinos. Note also that the case 
Tdec(^-dm) < Tdec(,y-e) might havc a significant impact on 
nucleosynthesis, so one has to treat it carefully. 
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3 WIMP DAMPING MASS 

To get quantitative estimates of the damping lengths, we 
need to specify the DM characteristics. In particular, we 
need to determine when DM decouples from neutrinos and 
electrons {i.e. taec(dm-i.) and taec(d,7i-e))- 

Similarly to Taec(dm~v), the decoupling temperature of 
DM from electrons, if mdm > T^ec^dm-e) > "1^, is given 

by: T, 



dec{dm — e) 



— — - — ' cm (where aa„ 



is expected to be temperature dependent with n = 2). How- 
ever, Tjec(dm-e) < "^e is Unlikely as this would require very 
large values of adm-e to compensate the fact that elec- 
trons become non relativistic and drastically annihilate be- 
low their mass threshold. 

More details are given in the next section but, gener- 
ically, the elastic scattering cross section of fermionic or 
bosonic DM candidates with neutrinos or electrons can be 
written as: 



Cdn 



where the DM couplings to neutrinos or electrons are "hid- 
den" in the proportional factor and with ms^, mse the 
masses of the exchanged particles (they are bosonic if DM 
is a fermion and fermionic if DM is a scalar). Note that 
mse (if charged) is necessarily larger than a few 100 GeV 
Eis no charged particle has been detected in past accelerator 
experiments, while msv can be below 100 GeV, provided 
one makes sure that these extra neutral particles are consis- 
tent with existing limits (they should be for example heavier 
than ~ 45 GeV if significantly coupled to the Z to avoid an 
anomaly in the Z decay width). 

When msv and mse > mdm, we find typically for a 
fermionic or bosonic DM: 

, in-54 f X"" 5-1 

beK.™-.) oc 10 (Y^^^j cm s 

, in-54 ( ^Se \-'^ 3-1 

6e,(.™-e) OC 10 [y^^^) cm S , 

leading to 



Td. 



'.ec(dm, — h') 



> MeV ( ) 
V 100 GeV / 



MeV/ ' 



Tdec(dm-v) and Tdec(dm-e) should actually differ by a fac- 
tor proportional to a certain power of the couplings but we 
do not expect Tdec(^dm-e) to be <C IMeV. Also these tem- 
peratures are expected to be slightly different depending on 
whether DM is bosonic or fermionic. 

There is an exception to these formulae, however, in 
the case of light DM particles (with a mass below 100 
MeV) which are required to get an annihilation cross-section 
dominated by a ter m in to satisfy gamma ray fiuxes 
iBoehm et al. 2002a^ : there may be other ways out but this 
is the simplest one. This implies, as explained in the subsec- 
tion 3.2, that &ei(dm-i/,e) depcuds rather on mdm instead of 
msv,e. 

In case of a mass degeneracy msu or mse ~ mdm, these 



decoupling temperatures become smaller, potentially mak- 
ing the mixed damping relevant. However, even if there ex- 
ists a mass degeneracy large enough so that one would a 
priori expect Tdec(dm-e) to fall below 1 MeV, DM is likely 
to decouple from electrons at T > (due to the drastic 
change in the electron number density after e~ annihila- 
tions) . 

The fact that b^K^dm-e) and b^K^dm-v) appear very close 
indicates that the DM thermal decoupling may be given by 
dm — V interactions instead of dm — e. It would be wrong, or 
at least dangerous, to focus only on dm-e interactions as the 
dm — V interactions determine whether the mixed damping 
regime exists or not. 

Before computing the damping lengths of a typical 
fermionic candidate (say neutralinos) and a scalar candi- 
date, we just mention that, in a situation where there is 
no mixed damping, the largest damping appears to be the 
free-streaming. The latter is indeed expected to erase all 
primordial fiuctuations with a mass below: 



M/s cx Mq c 



mdr 



-15/8 



ms 



'3/2 



MeV ) V 100 GeV / 
where c denotes the couplings (at a certain power). 

3.1 Neutralinos 

The expected ma ss for neutralino DM lies between 
[0(GeV),0(TeV)] llBoehm et al. 200(i) . The lower limit 
set by the LEP experiment however is ~ 37 GeV 
( ALEPH collaboration 2001) but this has been obtained by 
assuming a gaugino mass unification. Also within the Mini- 
mal Supersymmetric Standard Model (MSSM), the upper 
limit set by relic density calculat ions based on s tau co- 
annihilations is about 0(400) GeV jEllis et al. 200(i) . Thus, 
one may find interesting to look at the damping mass of pri- 
mordial fiuctuations associated with neutralinos in the mass 
range [0(40) GeV, 0(400) GeV]. (Looking at lighter neu- 
tralinos would actually be very interesting too, but if they 
are lighter than 10 GeV, one needs to make sure that their 
annihilation cross-section is suppre ssed by the square of the 
DM velocity to satisfy radio fiuxes iBoehm et al. (2002ah . a 
condition which can perhaps be achieved by introducing, for 
example, a new light gauge boson that is not included in the 
Standard Supersymmetric Model.) 

Let us therefore focus on neutralinos heavier than 40 
GeV. The damping induced by neutrinos is basically given 
by the interactions Xo ^ Xo through the exchange of 
sneutrinos. The latter has been searched in LEP experi- 
ments, notably from single photon events, and a limit of 
~ 84 GeV has now been set on its mass (under the assump- 
tions of gaugino and sfermion mass unificatio n as well as no 
sfermion mixing J ALEPH collaboration 2002M . 



The total elastic scattering cross section is given by: 



IGtt (m| — m'l)'^ ' 



(14) 



with ci the coupling between sneutrino, neutrino and neu- 
tralino. There are two cases for which it is interesting to 
compute the damping scale, namely: 
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• a heavy or light neutrahno much lighter than sneutri- 
nos, {mc, > m-^), 

• a heavy or light neutralino degenerate with sneutrinos 
(rac, — m^) so as to benefit from a resonance that would 
enhance the cross section and therefore the damping effect. 

Let us first consider a case where there is no mass de- 
generacy. One gets 

-4 

^el{dm — ly) 



100 GeV. 

which appears to be much smaller than the value required to 
have mixed damping. So free-streaming is the most relevant 
effect, erasing up to: 

,f 3/2 /'mdmX-i^/s / mc. \-3/2 



M.. 8 10- -10-^ 



-3/2 



Mp 



40 - 400 GeV. This is to be com- 
pared with 

" V0.3/ VlOO GeV 
for = 40 — 400 GeV. All these values are actu ally far 
from the calculation done by jHofmann et al. 200 ill . 

In case of a mass degeneracy, on the other hand, the 
Xo — I' elastic cross section, 

(-m? 4- ml 2rmx) 

can become temperature independent provided — < 
2T m^. Let us define rap as mp = m-^{l + 5). So 

• for T > (2 5 + 5'^) m^/2, the cross section is temperature 
independent: 



(15) 



4 m2 ' 



• while for T < {25 -\- 5'^) the cross section is tem- 

perature dependent: 



(16) 



For example, one expects (Jx-v to be constant at tempera- 
tures above 10.5 GeV, if one assumes m-x ~ 100 GeV and 
5 = 0.1. 

Thus, unless the degeneracy appears to be adequate so 
that DM can decouple at T < 1 MeV, the damping gener- 
ated by the dm-i/ collisions should eventually be computed in 
two steps. First using formula [T51 UP toT^ = {2S + S^) m;^/2 
and then using eq^Hup to Taec(dm-v) (where Taec(dm-u) will 
be estimated by use of ea ll6ll . One nevertheless expects the 
damping to be dominated by the late times so using ea ll6l 
only should finally give a good estimate. 

One can compute the neutralino thermal decoupling 
temperature in ordinary situations from ea ll4l using eq|S| 
and eq|5| One finds 

1/4 



dec{dm~^) ~ 5.3 MeV ( — 



nidv 



100 GeV J V MeV ) 



for T < (2(5 + S^)mx/2 (i.e. 75 MeV for ma^ = 40 GeV 
and 133 MeV for mdm = 400 GeV). A similar temperature 
is expected for Tdec{dm-e) as the electrons are relativistic 
above rrie and behave like neutrinos. On the other hand, if 
the degeneracy was about 10~^ for mdm = 100 GeV (i.e. 
5 = 10~^, mc = 100.001), then only one step would be 
necessary to compute the thermal decoupling epoch as the 
cross section remains constant up to T < MeV. As a result, 
Tdec{dm,-v) could be less than ~ 1 MeV for m^ ~ 100 GeV, 
which would further considerably increase the damping mass 
(at least up to 1 Mq), although such a degeneracy cannot 
be taken seriously. 



3.2 Light DM 

Although the most popular DM candidate certainly appears 
to be a "heavy" neutralino, the possibility of having light 
DM candidates may still be a viable solution. In partic- 
ular, if not coupled to the Z, light scalar particles (with 
O(MeV) < mdm < 0(10 GeV)) may be solution to the DM 
issue, provided they are coupled to a light gauge boson U 
(with a mass mu > mdm) or to heavy fermionic particles 
F (although in this case, the gamma ray, radio flux as well 
as <7 — 2 constraints might impose that DM be made of non 
self-conjugate particles and require to introduce a set of new 
particles in order to kill the too large F contribution to the 
muon and electron g ~ 2). 

When coupled to fermionic particles F, light scalars are 
not expected to experience a large damping effect because of 
the mass difference between DM and the F particles (sup- 
posed to be very heavy from past accelerator experiments). 
There is an exception, however, if DM mainly annihilates 
through a very light particle F° (with a mass mpa ~ mdm) 
which is not significantly coupled to the Z (and which would 
have escaped accelerator limits). 

On the other hand, if DM is coupled to a new and light 
gauge boson U , one expects damping effects to be signifi- 
cantly enhanced. The elastic scattering cross section Odm-v 
(associated with the exchanged of a [/ boson through a t- 
channel) is given by: 



"eZ(dm — ly) 



10" 



S^2 f2 
JU, 



mu 
MeV 



In iBoehm fc Favet f2003t) . it is shown that - to satisfy 
both relic density and g — 2 constraints - the product 
of the couplings Cu and fu^ (which correspond to the 
U boson couplings to DM and to ordinary particles re- 
spectively) must satisfy the relationship Cu fui = (3 — 
12) 10-«(^)-i (^)^ so 



'-'eZ (dm — f) 



[9 - 144] 10" 



md„ 
MeV, 

I.e. same value for U 



(where we assume "universal" fui 
e — e and U ~ /j. — jl for instance). This implies a decoupling 
scale-factor: 



Udecidm-,.) ~ (5 - 10) 10 



mdr, 



Mev) 



-3/4 



which is independent of mu and maximal for small values 
of the DM mass. As an example, the decoupling epoch for 
mdm ~ 10 MeV corresponds to a£jet,(tim-i/) ~ (0.8—1.5) 10~^ 
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which implies a mixed damping mass of about 

M„d ~ 21 - 172 Mq (for lOMeV). 

This damping mass obviously gets smaUer for larger DM 
mass and it is worth to mention that the cross section 
for mdm = 10 MeV is close to the upper limit that 
arises from the condition of free-streaming neutrinos, namely 

foeK.™-.) < 1.8610-^« x-^e^-- cm^s-i 

(so vridm = 10 MeV is close to the upper limit for &e;{dm-i/) ~ 
144 10"*^ (^f cm'^s"^). Note that one would naively 
expect a bigger effect in the case of light DM due to the 
smallness of mdm but the latter is actually compensated by 
the smallness of the product Cu fui ■ 

With the expression adec(dm-v) obtained above (valid 
only if DM particles are coupled to a [/ boson), one can 
readily see that the mixed damping regime is at work only 
for a DM mass in the range mdm £ [~ 5 , (500 - 700)] MeV 
(heavier DM would decouple before 1 MeV, lighter particles 
would give rise to the alternative scenario where the neutrino 
decoupling occurs at Tdec(v) < MeV). 

In the previous examples, however, we assumed that 
the number density of DM particles was equal to the num- 
ber density of anti DM particles. If one relaxes this assump- 
tion (assuming non self-conjugate DM), then Cu can take 
much larger values (but likely to be G [10"'', 0(1)]). As 
a result, the decoupling scale-factor which is proportional 
to b^(*dm-v) .fu,y^'^ can become at most~ V 10^ 

larger than our previous estimate, so that the damping 
mass that one expects in such situations could possibly 
reach ~ {VWf ~ 3 10* times the mass Mmd estimated 
without the number density asymmetry. But a change in 
Cu is, in fact, particularly relevant for DM particles with 
mdm ^ 0(10) MeV, for example, as the upper limit on the 
cross section (set by the condition of free-streaming neutri- 
nos) may remain possible to satisfy. 

Such particles, however, would not have any residual 
annihilations so they would not produce gamma rays but on 
the other hand, they would have an oscillating linear matter 
power spectrum largely damped below 0(10'' Mq), which 
finally provides a signature. 

One can also consider larger values of fui by relaxing 
the assumption of universal couplings. Indeed, the values 
mentioned before for Cu fui actually supposed that the cou- 
pling dm — u — u is of the same order of magnitude as the 
couplings dm — jj. — jj, or dm — e — e (constrained through the 
U contribution to the muon and electron anomalous mag- 
netic moments). Writing fu^ = xfu^ (with a:: > 1) {fu^ being 
the relevant coupling for low DM mass), one finds that the 
damping mass is increased by a factor x^^'^ and a(t) by a 
factor x^^'^ (still provided one makes sure that the condition 
of free-streaming neutrinos is satisfied). 

Such situations may finally allow for DM particles heav- 
ier than the limit 500 — 700 MeV to experience a non- 
neghgible mixed damping (while they were not supposed to 
with smaller values of the couplings) . Indeed, an increase in 
the couplings simultaneously allows for ea ll3l to be satisfied 
and for Mmd to be enhanced. 

Note that the mixed damping mass is not expected to 



become larger in case of a degeneracy between mu and mdm, 
because the interaction at the origin of the damping effect 
proceeds through a t-channel. On the other hand, such a 
degeneracy would affect the decay modes of the U boson 
which would then decay into ordinary particles only {i.e. 
into e^e~ and vi' instead of dm dm*, although one has to 
make sure this is not excluded by nuclear experiments). 

The previous examples provide a theoretical framework 
to investigate the effect of a "coUisional" cut-off in the linear 
matter power spectrum at scales ~ 10'^ Mq. This is actu- 
ally too small compare to our resolution but the question 
of whether a cut-off in the linear matter power spectrum 
is still present in the P(fc)ni can be partially answered by 
investigating the effect of a cut-off at 10^ — 10^ Mq. 

The main feature, say oscillations, in the light DM lin- 
ear matter pow er spectra should be somewhat similar to 
that obtained in lBoehm et al. f2002d) (although the damp- 
ing mechanisms are different). However, to extend our point 
to other class of DM candidates, we decide to perform simu- 
lations of coUisionless WDM and see whether a cut-off in the 
linear matter power spectrum at large scale (lO'' — 10^ Mq) 
also appear in the P{k)ni and at which redshift. If such 
a cut-off appears difficult to detect, then non conventional 
candidates as proposed in this paper will be very difficult 
to exclude from the measurement of the non linear matter 
power spectrum only. 



4 REGENERATING SMALL SCALE POWER 

Estimations of the matter power spectrum on scales I > 200 
kpc are now routinely derived from weak and strong 
lensing techniques at low redshift, or indirectly from the cor- 
relation of spectra of the Ly-a forest at 2; ~ 3. Fair agree- 
ment of the observations with the slope of the non-linear 
power spectrum expected in the CDM paradigm is com- 
monly regarded as a confirmation of the assumed "cold" na- 
ture of dark matter. This is a biased conclusion, as the n-ni = 
— 1.4 slope of P(fc)ni at these small scales (n = d In P/A In k) 
and at this low redshift range is known to behave as an "at- 
tractor" to the evolution of w , for a variety of DM models 
iScoccimarro &: Frieman 199dl . 

Using numerical simulations of structure formation with 
initial power spectrum restricted to two separated scale 
ranges, Bagla fc Padma nabhan ( 1997.) demonstrate that the 
non-linear growth of the small scale modes is mostly driven 
by the cascade of power due to the coupling with larger scale 
modes. As a result, DM models with very different P{k)\ 
on small scales but similar large scale power are expected 
to yield small-scale P{k)ni of the same shape. Focusing on 
measurements of the Ly-a forest, IZaldarriaea et ah 
make the same point at a somewhat more observational de- 
gree, and clearly show how little information is obtained on 
the initial shape of P{k)i at fc 10 /iMpc~^ by inverting the 
flux power spectrum, as long as one starts with a sufficiently 
large prior in ni. 

Physically, the IDM and WDM models considered here 
are as well motivated as CDM (but for a little more com- 
plexity), and they are a priori as probable as the latter. It 
is therefore necessary to verify if, when and at which scale 
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their specific matter power spectrum becomes sufficiently 
close to that of CDM so that current observations are un- 
able to separate between the two. 

To assess this, we simulate the gravitational formation 
of structure using 128"^ coUisionless particles in a small, 5 

Mpc side comoving box, with the cosmologic al param- 
eters measured by WMAP iSoergel et al. 2003^ ■ starting 
from z=100. The mean inter-particle distance is 39 kpc 
and our force resolution is one tenth this value. We evolve 
the particle distribut ion with the public version of GADGET 
JSprineel et al. 200 ill . We consider a CDM and four WDM 
models with warmon masses 0.6, 1.1, 2 and 3.5 keV: the 
associated free-streaming lengths correspond to Lagrangian 
masses 10®, 10*, 10^ and 10** h'^ Mq respectively. The 
CDM and 0.6 keV WDM models bracket the series of mod- 
els discussed in this work (including IDM), while the other 
3 samples are intermediate cases. We keep the same phases 
for all the modes in the five simulations. In the following 
discussion we will not tackle the issue of cosmic variance. 
We normalize the modes 5k inside the box so that they cor- 
respond to (Tg — 0.9. In doing so, we neglect the effect of 
primordial fluctuations with comoving scales greater than 
the box size: the power effectively realized in the box is 3 
percent of the total cosmological power. 

To estimate the dispersion expected if we were to in- 
clude all the modes up to the scale of cosmic homogeneity, 
we have performed two series of similar simulations which 
we call test simulations. Checking this effect is necessary 
because of the combination of our limited dynamic range in 
scale and of our choice of a small simulation box. 

The first series of test simulations consists of a CDM 
and a 0.6 keV WDM model with the same phases as above 
in a flat cosmology with the same and as and to but 
with different flo: f2o=0.15 and 0.45, corresponding to a 
^>5h~^ Mpc ~ ±""5/3 background and r2o=0.03 and 0.57, 
corresponding to a 5^^^-i jyj-p^, ~ ±o"5 background. In these 
simulations, h varies accordingly to ensure the same age to 
as the cosmological model preferred by the WMAP results. 
Recall that we take as = 0.9 and we define as as the RMS 
density fluctuations in spheres of radius 5 Mpc. This 
series of simulation checks the effect of large scale bias in a 
peak-background split approach. 

While we flnd the regeneration of small scale power for 
our 0.6 keV WDM model to be very similar in the f2o=0.45 
and ^0=0.57 simulations to the result of the reference mean 
density simulation, the low density simulations (especially 
no=0.07) show that small scale regeneration occurs signifl- 
cantly earlier when one evolves in an underdense, 5 Mpc 
region of the universe compared to a mean or high density 
region. This is expected, because for the same absolute nor- 
malisation of the power spectrum, the growth of modes oc- 
curs earlier in low density flat cosmologies. As a result, the 
cascade of power from the larger modes that we describe 
in detail below also occur earlier in low density cosmologies 
and the 0.6 keV WDM power spectrum catches up earlier on 
with the CDM power spectrum. However, this phenomenon 
will only strengthen our conclusions. 

The second series is the same CDM and 0.6 keV WDM 
model as in our main series, changing only the box size from 
5 to 20 Mpc (and the mean inter-particle separation to 



156 kpc). This second series verifles that enough large- 
scale power is already realized inside the 5 Mpc box 
for the same cascade of power to occur as would happen 
in a larger box. We have found that the neglected larger 
scale modes have little impact on our conclusions: the dis- 
persion due to fluctuations in matter density is smaller than 
the differences we discuss in the power spectra of different 
DM models, and the regeneration of small-scale power is 
the same in the 5 and 20 Mpc boxes at the overlap- 
ping wavelengths. For completeness sake we point out that 
we did flnd reduced small-scale regeneration in a 1 Mpc 
box simulation because of lack of large scale power in this 
case. This allows us to conclude that 5 Mpc is larger 
than, but close to the minimum box size necessary for our 
conclusions to be meaningful. 

Figure gives the matter power spectra P(fc)ni of our 
five simulations, from z=30 to z=0. The spectra are com- 
puted up to /cNy = 2 7r64/5 ~ 80 /iMpc~^ using a de- 
convolved TSC scheme. To facilitate the comparison, they 
have been divided by the normalized linear growth factor 

The qualitative evolution of the pow er spectrum is simi- 
lar to that found bv lKnebe et al. ^2002^ in their simulations 
of WDM models with warmon masses of 0.5 and 1 keV. 
At z=30, shortly after starting the simulations, the shape 
of the measured spectrum is in all cases close to the lin- 
ear power spectrum of the initial conditions. By z=10, non- 
linear effects have already significantly altered the shape of 
all power spectra on the smallest scales > 10 /iMpc~^. At 
z=7, WDM power spectra are already reasonably close to 
CDM, with the exception of the 0.6 keV simulation, where 
the power spectrum is still steeper and keeps the signature 
of the 10® Mq cut-off. At z<2 however, the power spec- 
trum of the 0.6 keV WDM simulation is within a factor of 2.5 
of that of the CDM down to the smallest scales probed. At 
z=l, the curves can be considered the same within 10 per- 
cent. These plots confirm that the late P(fe)ni is similar down 
to scales as small as 80 kpc at z<2 for models with so 
different input power spectra such as CDM and the 0.6 keV 
WDM scenario. We find the slope of the evolved power spec- 
trum to tilt from n ~ —1 for A; < 10 hMpc~^ to n ~ —1.6 
for fc > 10 /iMpc"^. Furthermore, at a fixed comoving scale, 
the power spectrum of WDM models with the larger warmon 
mass (and smaller free-streaming length) catches up earlier 
with the CDM power spectrum than do models with larger 
free streaming lengths. The intuitive guess that at a fixed 
WDM model, power on large scales matches earlier CDM 
power than power on a smaller scale is confirmed. This is 
of course the signature of the exponential cut-off present in 
the linear power spectrum. 

These monotonic behaviours are more clearly seen in 
Figure |21 There, we first divide our wavelength coverage 
(fe € [0.1 1.9] /iMpc~^) into 6 consecutive bins, correspond- 
ing to the 6 panels shown. In each bin we estimate the aver- 
age of the power spectrum (P(fe)ni) for the CDM simulation 
and the four WDM simulations at each redshift. We then 
divide the time evolution of the averaged power spectrum 
{P{k)ni{z)) of each WDM simulation by the time evolution 
of the averaged CDM power spectrum. We only discuss the 
upper left and lower right panels, the other 4 being inter- 
mediate cases shown here for reference. On the upper left 
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panel (bin associated to the largest scales), the ratio of the 
average power measured in the WDM models to the aver- 
age power in the CDM model is within 10 percent of unity 
from z=30 down to z=0 for all warmons but the lightest, i.e. 
the 0.6 keV particle. In this last case, the signature of the 
initial exponential cut-off is already visible on large scales, 
with a factor 2 less power than CDM. This deficit gradually 
disappears over 10>z>2. The lower right panel (bin associ- 
ated to the smallest scales that we probe) clearly separates 
our four WDM models. The average power of the 3.5 keV 
WDM model is the first to catch up with CDM, at z~12. 
The increase in the ratio of 2 and 1 keV WDM power to 
CDM power is abrupt before z~10, slows down in 10>z>2 
and becomes comparable to CDM at low redshift. The ratio 
of the 0.6 keV WDM to CDM power increases less steeply 
than the 1 and 2 keV power initially, seems to slow down at 
z<l, but has not yet exactly matched the CDM power at 
z=0: there is a remaining 10 percent offset (but probably in- 
distinguishable with current observations). However, in all 
these cases, and also including the four remaining panels, 
the "monotonic" feature of the process of small scale regen- 
eration is preserved. Figure|5|contains the same information 
as Figure it only has more sampling points in redshift. 
We can extrapolate that the power of the 0.6 keV model 
may not have matched the power of the CDM model at 
z=0, if one would look at scales smaller than logj^Q k = 1.9. 
As we want to simply illustrate the late-time degeneracy in 
the power spectrum between a 0.6 keV WDM model and a 
CDM model at scales ~ 100 hT^ kpc, we leave this issue for 
further work. 

Regarding the regeneration of small-scale power we find 
that observations (e.g. lensing) probing the matter distribu- 
tion at z<l at scales as small as 100 h"^ kpc will be unable 
to detect the signature of the initial exponential cut-off of 
WDM models with warmon mass as small as 0.6 keV. For 
this warmon mass, the signature of a cut-off will only be 
clearly visible on lenses at 7.<^ 4, not even considering sys- 
tematic errors nor cosmic variance. All WDM models with 
larger warmon mass will be increasingly harder to single out. 
Alternatively, one may need to probe scales smaller than 100 

kpc. The converse, and our main point, is that find- 
ing at the above scales a z=0 power spectrum in agreement 
with that expected from the non-linear evolution of an initial 
CDM model cannot rule out physically motivated alterna- 
tives such as WDM or a light DM scenario presented in this 
paper. 

We emphasize that this is true even if we take into ac- 
count the possible effect of the modes with wavelength larger 
than our simulation box: while evolving in a 1 cr overdense 
region (defined with respect to our neglected modes) does 
not change the regeneration of small scale modes compared 
to what happens in a mean density region, simulating a 1 a 
underdense region we find small scale regeneration to occur 
even earlier. 

To conclude this section, we note that the mass func- 
tion of collapsed objects is also a sensitive probe at high 
z of a possible exponential cut-off in the primordial power 
spectrum. It has the advantage over power spectrum mea- 
surements that large variations in the abundance of haloes at 
the epoch of the formation of the first stars or later at reion- 
ization can be rather tightly constrained even with current 



observations, at a fixed physical prescription for the evolu- 
tion of the baryonic component. The drawback resides in 
first, the current large freedom in physical models of reion- 
ization or of the formation of the first stars and second, 
in that the mass function is also very sensitive to a possi- 
ble small-scale primordial non-gaussianity. While we plan to 
address these issues in more detail in the future. Figure 
gives the z=0 (left panel) and z=10 (right panel) mass func- 
tion of collapsed dark matter haloes in our simulations. Note 
the mass range probed (5 x 10^ to 10^^ Mq) At z=10, 
the deficit of haloes of all masses is evident in the 0.6 keV 
WDM model compared to CDM. It is also apparent in the 
simulated WDM models with higher warmon mass, although 
mainly at the low-mass end of the curves. The z=10 abun- 
dance of the largest, 10^" Mq mass haloes is similar in 
the CDM and the 1.1, 2, 3.5 keV WDM models. At z=0 
the abundances of haloes with Mtot ^ 10^^ Mq match 
in all 5 models. In the range 10«-5 < Mtot < 10"/i-' Mq, 
the mass function of the 0.6, 1.1 and 2 keV models shows a 
deficit of haloes compared to CDM, with the 0.6 keV model 
showing the lowest abundance. At Aftot < 10^'^ Mq, 
halo abundances tend to coincide again. 



5 CONCLUSION 

In this paper we compare two distinct DM candidates. One 
is the most popular candidate, say the lightest neutralino, 
while the other one is light DM particles. 

Both these candidates are expected to have weak inter- 
actions. They both annihilate to get the proper relic density 
but, in the case of neutralinos, coUisional damping effects 
are completely negligible while, in the case of light DM, the 
coUisional damping (and more specifically the mixed damp- 
ing) yields a cut-off in the linear matter power spectrum 
which can realistically reach 10^ Mq . 

Because the existence of this cut-off seems a promising 
way to discriminate among these two extreme DM scenarios, 
we perform a series of numerical simulations to obtain the 
corresponding non-linear matter power spectra. To make the 
case even more obvious, in fact, we simulate scenarios with 
a more drastic cut-off than what is expected in our model. 
More realistic simulations of light DM scenarios will be pub- 
lished in a forthcoming paper but we do expect similar con- 
clusions. This enables us to quantify the time evolution of 
P{k)ni as a function of scale, and also to precisely assess to 
which extent the late time (z < 1) small-scale power spec- 
trum can be used as an efficient tools to distinguish between 
standard (coUisionless) CDM and alternative scenarios. 

Since we are limited by the size of our simulations, we 
simulate a cut-off at 10'^, lO'', 10® and lO" h'^ Mq. We 
show that the corresponding P(k)ni is similar to the non- 
linear matter power spectrum of CDM particles for z < 1, 
so weak lensing measurements except on the smallest scales 
(100 kpc/h or less for lenses located at z < 1) are not able 
to discriminate between usual CDM particles and more ex- 
otic candidates. Note here that current measurements of the 
cosmic shear probe the shape of the matter power spectrum 
down to 200 kpc, while planned upcoming surveys will 
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Figure 1. Snapshots of the dark matter power spectra measured in simulations of a series of DM models. CDM and 0.6, 1.1, 2 and 
3.5 kcV heavy warmon WDM corresponding to respectively 10®, 10*, lO'' and 10® Mq Lagrangian masses in the free-streaming 
length are shown. The power spectra have been divided by the linear growth factor to facilitate the comparison between redshifts. Note 
the exponential cut-off at z= 30 in the spectrum of WDM models and the similarity of the spectra of all five models over the simulated 
scales at z<2. 



go down to 70 ft ^ kpc and start to probe a critical region 
for a 10^ h-^ Mq cut-off. 

On the other hand, there exist other sensitive signa- 
tures of the nature of the dark matter, such as structural 
parameters, clustering of virialized dark matter haloes and 
the reionization epoch. However these are more indirect con- 
straints and involve poorly understood interactions between 
baryons and DM, which considerably weakens their predic- 
tive/discriminative power. 
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Figure 2. Evolution with redshift of the ratio of (1) the mean of the power spectrum of each of the four WDM models (0.6, 1.1, 2 
and 3.5 keV particles corresponding to respectively 10^, 10®, 10^ and 10^ Mq Lagrangian masses) to (2) the mean of the power 
spectrum of CDM model, in six consecutive bins covering the whole wavelength range probed by the simulations. Note how the ratio 
converges to unity at late times for all WDM models in each of the 6 bins. The marginal exception is the evolution of the 0.6 keV WDM 
model which shows 10 percent less average power than CDM on the smallest scales bin at z=0. The power spectrum of this model may 
strongly depart from CDM on scales logiQ{k) ~ 2 even at z=0. 
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Figure 3. Early (z=10) and late-time {z=0) dark matter halo mass functions of the CDM model compared to those of the WDM 
models. Note the deficit by a factor of ~ 40 of 0.6 keV WDM haloes at z=10 compared to the abundance of CDM haloes. 



